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We propose a novel approach in the study of transport phenomena in dense systems or systems 
with long range interactions where multiple particle interactions must be taken into consideration. 
Within Boltzmann's kinetic formalism, we study the influence of other interacting particles in terms 
of a random distortion of energy and momentum conservation occurring when multi-particle inter- 
actions are considered as binary collisions. Energy and momentum conservation still holds exactly 
but not in each model binary collision. We show how this new system differs from the Boltzmann 
system and we note that our approach naturally explains the emergence of Tsallis-like equilibrium 
statistics in physically relevant systems in terms of the long since neglected physics of interacting 
and dense systems. 
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The Boltzmann equation (3] for a spatially independent 
system of particles describes the evolution of the one par- 
ticle momentum p distribution f(pt) of a "foreground" 
particle with mass m subject to interactions (here as- 
sumed to be number conserving) with N "background" 
particles j = 1, . . . , N with masses mi, W2, . . . , mjv- In 
classical, non-relativistic regimes, the Boltzmann equa- 
tion takes the form 



[W(p',p)f(p', t) -W(p,p')f(p, t)] d 3 p', (1) 



dfjpt) 
dt 



where we have suppressed a possible time dependence in 
the transition rate W(p,p'), the rate per unit time of the 
foreground particle making a momentum transition from 
p to p' due to interactions with the background particles. 

In applications, the Boltzmann equation is restricted to 
sufficiently rarefied systems with short range interactions 
so that only two particle interactions (two body or binary 
collisions) are incorporated into the transition rate |^|. 
For these systems, the transition rate is given by || 

W% p ')=£ [ / d% f d%' S 3 (p + q- - p' -q/) 

x5(E p + E 9j . - E„/ - E q >) (Hi h (#,*)] . ( 2 ) 

where the energies are given by E p = \p \ 2 /2m and E 9j = 
\q~j\ /2mj, bj(q~j,t) is the momentum qj distribution for 
background particle j, and (crv)j is the differential cross 
section for the scattering of the foreground particle with 
a background particle j. 

No satisfactory generalization of the transition rates 
has been given to describe transport phenomena in dense 
systems or systems with long range interactions, that is, 
in systems where multiple particle interactions must be 
taken into consideration The natural generalization 
of the transition rates is to assume that the transition 
rate is proportional to the sum over collections of back- 
ground particles of 1) the expectation of a collision be- 



tween the foreground particle and a collection of back- 
ground particles and 2) the expectation of the conser- 
vation of energy and momentum constraint. However, 
there is no clear way of systematically generalizing the 
Boltzmann equation by adding successive multi-particle- 
interaction correction terms (which are generally diver- 
gent). Attempts to generalize the Boltzmann equation 
for even "moderately" dense systems have largely been 
abandoned. 

In dense systems or systems with long range interac- 
tions, for the transition p — > p ' of the foreground particle 
to be physically possible in an interaction, we need the 
initial and final momenta and energies to satisfy 
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by the conservation of momentum and energy. Singling 
out the dominant role that a single background particle 
j plays in making the transition physically possible, we 
write the above constraints as 



P + Qj = P +Qj +e 
E p + E 9j = E p ,(2 7 -1) 



E„ 



(5) 
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The latter form of the constraints appears arbitrary. 
However, when combined with Eqs. (|3j) and (^), it serves 
as a defining relation for the variables e and 7. The above 
form is chosen here as such for mathematical convenience. 

Thus, in dense systems or systems with long range in- 
teractions, for the transition to be "due" to an interaction 
with particle j, we need a foreground-background binary 
interaction as before, but with modified energy and mo- 
mentum constraints. The transition rates for dense sys- 
tems or systems with long range interactions are then 
given by simply replacing the <5-functions appearing in 
Eq. (0) according to: 
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S 3 (p + qj-p' -%') -» <(<5 3 (p + <Zj - f -<%' -e)))g 
5(E p + E Qj -E pl -E ql ) -> «*[E P +E a -E p ,(27-1)-E 9 j]», 



(7) 

(8) 



r 



So far, we have pushed all of the complexity of comput- 
ing multiple scattering transition rates into computing 
the distributions of the random vector e and the random 
variable 7. However, the form of the transition rate is 
now familiar and many aspects of the distributions are 
immediately clear from their definitions. As the density 
of a system and/or the range of interactions decrease, 
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for only two particle collisions are increasingly present. 
Clearly, our transition rates simplify to the usual binary 
transition rates for rarefied systems with short range 
interactions and thus, incorporate the successes of the 
Boltzmann equation in describing transport phenomena 
in these systems. 

However, as the density of a system and / or as the range 
of interactions increase, Var(7) and Var(e) increase. For 
such physically relevant systems, the success of our pro- 
posed generalization to the Boltzmann equation will be 
made evident by indicating how our generalization de- 
termines generalized equilibration dynamics which dif- 
fer from canonical equilibration dynamics and how our 
generalization very naturally reproduces experimentally 
measured equilibrium transport phenomena. 

In an interaction involving a foreground particle and 
a background particle j we can interpret the presence of 
other interacting background particles in dense systems 
or systems with long range interactions as providing a 
mechanism for carrying off or providing the additional 
momentum and/or energy necessary to make transitions 
physically possible. We can also interpret the above tran- 
sition rate as the usual transition rate in a system in 
which only binary collisions occur which do not neces- 
sarily conserve kinetic energy and momentum, as is the 
case for inelastic collisions. 

Not only do these two interpretations aid us in under- 
standing the nature of the 7 and e distributions, but we 
see that by extending our original definitions of e*and 7, 
the form of our transition rates is applicable to a very 
large class of systems: with elastic and/or inelastic inter- 
actions (here defined generally to incorporate all inter- 
actions which do not conserve energy and momentum), 
with low or high number densities, with short and/or long 
range interactions, and with any mechanisms which can 
carry off or provide additional momentum and energy to 
particles. 

Indeed, for general systems, the calculation of the e and 
7 distributions has become even more complicated. Be- 
fore proceeding, it is helpful to note that the expectation 
of the parameters of our generalization appearing in the 



delta functions of Boltzmann's formalism has the effect 
of widening or broadening the delta functions. Thus, the 
existence of our parameters with some fluctuation is guar- 
anteed in that no exact delta function exists in nature. 
As such, the importance of our proposed generalization 
to physically relevant systems depends on the extent to 
which systems exhibit fluctuations in the parameters 7 
and e. 

It is natural to expect that systems with large num- 
ber densities, long range interactions, inelastic interac- 
tions, and other such mechanisms (such as quantum me- 
chanical interactions) exhibit sufficient fluctuations to, 
and indeed as they are known to, require a generalized 
formalism. Thus, to verify the necessity of generaliz- 
ing Boltzmann's formalism to a system requires exten- 
sive theoretical calculations and/or precise experimental 
measurements of the 7 and e parameters of that system. 
One such theoretical calculation has already been given 
by H. Haug and C. Ell whose derivation of a semiclas- 
sical Boltzmann equation for Coulomb quantum kinetics 
in a dense electron gas suggests, in the limit of completed 
collisions, an asymptotic and approximate distribution 
of gamma given by a peaked distribution about 7 = 1: 

/ 7 ( 7 ) = 2 (r/2E;j / [( 7 - l) 2 + (r/2E;) 2 ] , where T is 
the sum of the collision dampening coefficients. 

In the absence of other such calculations or measure- 
ments of the distributions of 7 and e in various physical 
systems, we may still proceed in the verification of our 
theory since the precise form of our proposed generaliza- 
tion can be verified by comparing the consequences of our 
proposal with well established and currently unexplained 
phenomena. For instance, we can study transport phe- 
nomena within our formalism which result from certain 
classes of e and 7 distributions which we expect to be 
present in many physically relevant systems. 

In many systems, we expect e and 7 to be peaked dis- 
tributions about e = and 7 = 1 and approximately 
independent of the incoming and outgoing momenta and 
energies (so that we interpret the influence of other inter- 
acting particles, unaccounted interaction processes and 
mechanisms as producing a random noise distortion of 
the conservation of energy and momentum constraints). 
In particular, we might expect e to have a spherically 
symmetric distribution about e = 0. Moreover, since 
the family of Gamma distributions has a rich variety of 
shapes and can approximate many classes of distribu- 
tions, in many systems we expect the peaked distribution 
of 7 to be well approximated by a Gamma distribution 
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with parameters a and A, chosen so that the distribution 
is peaked about 7 = 1: 



/ 7 ( 7 ) = A(A 7 ) 



a-l -A7 



/T(«), 



(9) 



where the average and variance of 7 are Avg(7) = 
« 7 » 7 = a/A and Var( 7 ) = (( 7 2 )) 7 - «7» 7 = a/A 2 . 

The form of the distribution of 7 as a Gamma distribu- 
tion is chosen here for mathematical convenience and to 
analytically recover the precise form of Tsallis' proposed 
equilibrium distribution in the following example, but we 
emphasize that any peaked distribution of 7 about 7=1 
will yield a Tsallis- like distribution, i.e., a deviation from 
the Boltzmann equilibrium distribution. Numerical cal- 
culations with the distribution of gamma suggested by 
the work of Haug and Ell resulted in Tsallis-like equilib- 
rium distributions for a large range of T and the Boltz- 
mann equation was recovered as the width T — > 0. 

To analytically illustrate some immediate results of 
our proposed generalization, we consider a system of N 
classical background particles denoted by j = 1,...,N 
with mass m,-. Suppose further that the interactions cos 8 



are completely elastic and that the number density of 
the system and the range of the interactions are suffi- 
ciently small to guarantee that the momentum distribu- 
tion of the background particle j = 1, .. . , N approaches 
the Boltzmann equilibrium distribution with a common 
temperature T (measured in units of energy): b° q (qj) — 
Cj exp (— qj/2m,jT) where Cj is a normalization con- 
stant. After the N particles are sufficiently equilibrated, 
we inject a foreground particle with mass m 3> rrij for 
all background particles j — 1, ...,7V and an initial dis- 
tribution f(p,t = 0). Suppose that multi-particle col- 
lisions between the foreground and background parti- 
cles occur (for instance, as the diameter of the fore- 
ground particle is much larger than the diameter of the 
background particles or if the range of interaction is 
larger for a foreground-background interaction than for a 
background-background interaction) and/or that the in- 
teractions are inelastic. Thus, the distributions of 7 and 
e exhibit variation. Further, suppose (av)j is constant as 
in the case of hard sphere interactions. We also suppose 
that the background particles are sufficient in extent (i.e., 
the background is an ideal heat bath) so that we may 
assume that the background approximately remains in 
equilibrium during the evolution of the foreground par- 
ticle. Therefore, the time dependent distribution f(p,t) 
satisfies Eq. ([j]), where the transition rate is 

WW) = E/ d ^ J d % '((t 3 (P + % V q'j - e)»* 
j 

x<(«5[E p +E 5i -E p ,(27-l)-E^])) 7 

x (av) 3 b?{%). (10) 

Interchanging the order of integration and using the S 3 - 
function to integrate out qj ' (so that q.j ' = p+q—p' — e), 



we obtain 

wW)= 

P 12 \p — p' — £\ 2 W ~ p' — e\Qj cos 
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where we have introduced spherical coordinates 
{qj, 9j, 4>j) for the integration over qj. Note that we have 
chosen Qj to measure the angle between q and p — p' — e 
so that the integration over <j)j is trivial and yields a 
factor of 2tt. We can now integrate over cos9j using 
the remaining delta function. A required transformation 
introduces a factor (\p — p' — e \ qj)~ 1 mj and the single 
zero of the delta function uniquely determines cos 0j : 



p2 _|_ p/2 



P \p — p' — t\' 



2m 

The result is 



\P-P' - Alj' 
(12) 



in/):// , l^j^y-p^, ^ijj (13) 

where the integration over qj is taken over all positive qj 
which satisfy — 1 < cos 0j < 1. Solving this constraint 
implies that 



\p-p'-e\< 



Qj > ~ ~ A + I V'^l 
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to first order in rrij/m. The remaining integral over qj is 
trivial and one obtains 



E p + E p / 
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The final form of transition rate indicates that the dis- 
tribution of e determines the rate of equilibration, while 
the distribution of 7 determines the shape of the result- 
ing equilibrium. Now, if e is symmetrically distributed 
about e = (as we suggested above), we have that 



\p-p'-^ 
8m jT 

\p- p' - e\ 



8m jT 

\p-p' + A 



(16) 
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Multiplying Eq. (|T|) by ((exp(-jEjT)})^ interchang- 
ing p <-> p ' and making use of Eq. mm, we see that 

W(p<,p)(^e~^Jj =W(p,p')(^e~^Jl . 

V) 

Note that this is a detailed balance equation H , implying 



// ^ E P' 
Hip) = c((e T 



(18) 



is a time- independent solution of Eq. (Q), where C is a 
normalization constant. This is easily verified by sub- 
stituting f cq {p) in Eq. flU) for /(p,i) in Eq. (0) and by 
making use of Eq. (^). If the distribution of 7 is well 
approximated by a Gamma distribution with parameters 
a and A (as suggested above), we find that 



J </'W(l + ^A 



(19) 



Thus, interpreting the Tsallis 
eter as 



nonextensitivity param- 



<?t — 1 H — 

a 



(dm 



(20) 



and the inverse temperature of the foreground particle 
F = 1/T F as pp = f/3 = ((7)) 7 /3, we have that the 
equilibrium distribution of the foreground particle is a 
Tsallis equilibrium distribution with nonextensitivity pa- 
rameter gx, inverse foreground temperature f)p = 1/Tp, 
and normalization constant C: 



n(p)=C{l-P F (l-q T )V p ] 



(21) 



The Tsallis equilibrium distribution is often successfully 
used to model equilibrium distributions exhibiting power- 
tail behavior. Furthermore, Tsallis-like distributions are 
extensively measured in many physical systems, espe- 
cially in systems with long range interactions and in tur- 
bulent flows @. 

We note that our approach differs from Q in that our 
"natural" generalization does not up-front introduce ex- 
act Tsallis equilibration and, moreover, we retain the in- 
tuitive statistical form of Boltzmann's molecular chaos 
hypothesis, transport equation, and transition rates. Our 
proposal compliments the approach taken by C. Beck ||, 
which generates Tsallis equilibrium statistics by consider- 
ing fluctuations in the parameters of the Langevin equa- 
tion. We incorporate the approach taken by G. Wilk and 
Z. Wlodarczyk which generates Tsallis equilibrium 
statistics from fluctuations in background temperatures 
in the Boltzmann exponential factor, by incorporating 
these fluctuations into our distribution of 7. 



Concluding, we have proposed a generalization of 
Boltzmann equilibration dynamics to model transport 
phenomena in general, non-equilibrium systems. Our 
generalization resulted from our interpretation of the in- 
fluence of other interacting particles as a random dis- 
tortion of energy and momentum conservation occurring 
when multiple interactions are considered as binary col- 
lisions. As it turns out this also explains the appearance 
of Tsallis distribution in Fokker-Planck dynamics |pl| , 
where local energy conservation was not maintained. 

In outlook, we noted that our formalism is applicable 
to many different and large classes of systems: with elas- 
tic and/or inelastic interactions, with low or high num- 
ber densities, with short and/or long range interactions, 
and with any mechanisms which can carry off or provide 
additional momentum and energy to particles. We also 
indicated how the result of our proposed generalization 
can be determined in systems by studying classes of rea- 
sonable distributions for the proposed parameters of our 
generalization. Although we have only considered a sim- 
ple example where we could easily and explicitly deter- 
mine some immediate results, our approach is easily ex- 
tended to more general systems. We already have similar 
results for a relativistic, quantum mechanical foreground 
particle scattering in a thermal background of light rela- 
tivistic particles, and it is clear that the general situation 
is similar. 
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